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The basic building block for Lorentz invariant and ghost free massive gravity is the square root of 
the combination g~^r], where g~^ is the inverse of the physical metric and 77 is a reference metric. 
Since the square root of a matrix is not uniquely defined, it is possible to have physically inequivalent 
potentials corresponding to different branches. We show that around Minkowski background the 
only perturbatively well defined branch is the potential proposed by de Rham, Gabadadze and 
Tolley. On the other hand, if Lorentz symmetry is broken spontaneously, other potentials exist with 
a standard perturbative expansion. We show this explicitly building new Lorentz invariant, ghost- 
free massive gravity potentials for theories that in the background preserve rotational invariance, 
but break Lorentz boosts. 


Introduction. Much progress has been recently made 
in understanding massive gravity. Attempts to give a 
mass to the graviton date back to the work of Fierz and 
Pauli (FP) in 1939 [I|. They considered a mass term for a 
spin two field which is uniquely determined requiring the 
absence of ghost degrees of freedom (dof). The mass term 
breaks the diffeomorphism invariance of General Relativ¬ 
ity (GR), leading to a graviton with five degrees of free¬ 
dom instead of the two of GR. In 1972, Boulware and 
Deser (BD) found that the ghostly sixth mode, removed 
by the FP tuning, reappears at non-linear level or around 
non-trivial backgrounds 0 . It is now well known that it 
is possible to avoid the presence of the BD ghost by choos¬ 
ing a suitable potential Q . Such a ghost free theory can 
be constructed by adding to the Einstein-Hilbert action 
a potential V function of the symmetric polynomials of 
the square root of the matrix = g^°‘fau, where 
is a fixed reference metric Q. In the following we will 
consider the case = 77 ^^, i.e. the reference metric co¬ 
incides with the Minkowski metric. This theory, dubbed 
de Rham-Gabadadze-Tolley (dRGT) massive gravity, is 
Lorentz invariant (LI) in the unitary gaugeQ. More in 
general, the whole class of rotationally invariant massive 
gravity theories with five degrees of freedom can be con¬ 
structed ill by using the canonical analysis, which can 
be extended to any non-derivative modified gravity the¬ 
ory [|. 

In this work we present alternative consistent branches 
for LI massive gravity. We show how their perturba¬ 
tive expansion and construction is related to the spon¬ 
taneous breakdown of the Lorentz symmetry. We also 
briefly comment on some of their phenomenological con¬ 
sequences. 

Ghost free potentials. From the canonical analysis 
of a generic, non-linear massive deformation of GR in 


^ Invariance under coordinate transformation can be restored by 
introducing a suitable set of Stiickelberg fields Q . 


four dimensions Q 

S = J d^x^Mh {R -m^V)= Seh + Sy , (I) 

it follows that, in the unitary gauge, in order to have 
5 propagating dof, the potential V has to satisfy two 
conditions. In terms of the ADM variable 
they read 

RankI 1 = 3, V, + 2 i (2) 

where we define V = N V. Va is the derivative of V 
with respect to = Na = {N, W) whose components 
are the lapse and shifts; is the eigenvector associated 
to the null eigenvalue of Vab- Requiring that V has a 
residual Lorentz invariance besides rotational invariance 
forces the potential to be a function of the eigenvalues Xi 
of A = g~^r] i- We could consider a more general ref¬ 
erence metric in A, or even promote it to be dynamical 
in the context of bigravity. 

In the simple case of a 2d space-time, we can solve the 
partial differential equations in ([2]) under the assumption 
of Lorentz invariance and find that only two potentials 
are allowed Q 

V± = uqT®! ^\/Ai i \/A2^ = uoToiTr ^v^A . (3) 

With VX |± we denote the two possible different 
branches of the square root of X in two dimensions (mod¬ 
ulo an overall sign). The two potentials are therefore as¬ 
sociated with the two different branches of the matrix 
square rooi0. Extending this observation to four dimen- 


^ In general, an n x n matrix with n distinct non-zero eigenvalues 
has 2^ square roots, according to each of the possible choice in 
the sign of the square root of its eigenvalues. Assuming that a 
matrix M with eigenvalues Aj can be diagonalized by a matrix 
U, M = U we have that x/M = U yjU~^ where 

yjMif) = diag(di -v/AT, ... , di x/Xn)- We will assume that y/)^ 
are real and positive. 
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sions, one can check that a large class of potentials satisfy 
the conditions ([2]); 

^br = E ^ Ibr , (4) 

n\ 

n—O 


where 5'„ |br are the symmetric polynomials of \fX |br 
once a choice of branch for the square root is made 

5*0 = 1 , Si=Ti, S2 =tI -T2, 

S 3 = rf - 3 Ti r 2 + 2 T 3 , (5) 

with r„ = Tr (^y/X |br ^ The subscript |br implies that 
different options for the square root of X can be cho¬ 
sen, according to the ± sign in front of each of the 
square root of the eigenvalues. On the other hand, once a 
choice is made, the symmetric polynomials should be con¬ 
structed consistently and no mixing between the different 
branches is possible without violating the conditions ([2|) , 
i.e. without reintroducin^he 6*^ ghost modc0. 

The dRGT potential Q corresponds only to a sin¬ 
gle branch among all those in o, i.e. the branch 
where \JXy) = diag(-|-VAi, -I-VA2 , -bVA3 , -bVA4) = 
'JXd I dRGT. The other choices represent different poten¬ 
tials through we can realize non-linear ghost free theories 
of massive gravity. 

In the following we present two different approaches 
to study different branches for massive gravity. The first 
one is based on the construction of y/X by using ADM 
variables. Alternatively, we can construct y/X perturba- 
tively around a chosen background. 


quadratic equations 



AB-bB A = 


-I e>)' *' = 

f 0 e 


0 0 

0 ^(3x3) 


( 8 ) 


where 1C = = 7®-! — ff . We can focus for instance 

on two different branches for A parametrised as follow 


A - (I 

Vl - Q'^Q^ V -2' -2*2' 

where e? = 1 and 


0 0 \ 
0 V^lbrj 
( 9 ) 


Qf=eK^/^^ (/C- 1/2 lb,)*' e, (10) 

with The various branches of the square 

root of X are controlled by ca and by the branches of y/lC 
which is kept formal and what follows applies to any of 
its branches. Other choices of square roots are available 
besides this simple one, however they typically lead to a 
spontaneous breaking of rotational invariance. 

Notice that 0 is the most general transformation 
that trivializes the requirement det(VAB) = 0 (Monge- 
Ampere equation 0), which is the key property of po¬ 
tentials with 5 dof. Such a transformation has to be 
invertible due to the fact that the equations of motion 
already impose a constraint on the shifts and such re¬ 
quirement reads 


ADM approach. The ADM approach joins together 
the definition of y/lC given in 0 and the auxiliary vari- 
able used in Q to solve the first of ([2]). Following [T^, the 
square root can be expressed separating the dependence 
on the lapse: 

Ny/X = A + NM, (6) 

where A and B are independent from N. Squaring the 
above equation one gets a quadratic set of equations for A 
and B, leading to different possible branches. The choice 
made in [I^ corresponds to the dRGT potential; on the 
other hand different choices are possible. Let us show 
how different branches arise and their properties. It is 
useful to introduce a new shift variable such that, in 
the new variables, the action becomes linear in the lapse. 
The transformation is written implicitly as 

M = NC + Q\ (7) 

where Q* = Qf{"fjk, will be specified later. Solv¬ 
ing equation ([5]) for A and B translates in the following 


det(A^5/+55,Q*) ybO. (II) 

Once one finds a consistent background and perturba¬ 
tions around it are considered, (HU has to be satisfied 
order by order to get a standard expansion. 

Take now a flat background for g. Though both the 
reference metric and the background physical metric are 
Minkowski space, in the unitary gauge we allow the pos¬ 
sibility that they are not aligned, namely 

= diag(-c^I,l,I). (12) 

The misalignment is measured by c. Clearly, by a coordi¬ 
nate transformation one can transform (HU in 77^,., how¬ 
ever the reference metric will not be Minkowski anymore 
when c yf: I. Thus, even though V is Lorentz invariant, 
the background value of metric breaks “spontaneously” 
the Lorentz invariance SO(3, 1) while the rotational in¬ 
variance SO{3) is preserved. At the background level, for 
the ADM variables, we have N = c, the spatial metric 
7b = ^ij the shifts are zero, iV* = 0. 


^ Indeed one can check that, while are satisfied for any choice of 
sign among ±\/A7, the very same equations are not satisfied when 
the combination of two branches is considered. For instance, 
taking Vbri + ^ ^br 2 5 have 5 dof requires 6 = 0. 


• If is zero at the background level (5® = 0), to 
preserve rotational symmetry we need 

y/iC |br= VT |br= ei I3 , 


( 13 ) 






3 


with ef = 1. Hence 

VX = diag(eA c“\ ei l,ei l,ei 1). (14) 

Eq. (fTTl) gives 

c + caei ^ 0. (15) 

When c = 1, Lorentz invariance is not broken and ea = 

El- Thus '/X = EA VX IdRGT and this is just the dRGT 
branch up to an irrelevant overall sign. For c 7^ 1, Lorentz 
invariance is broken, then eq. (Ha is satisfied when ea = 
±Ei. In this case the choice of the negative sign is not 
the dRGT one and is non-linearly connected with the 
new branch: 

V Xd I New — diag(- a/aT, \/A, \/A 3 , x/A) , (16) 

this choice produces a genuine new potential. 

• If is not zero, computing the inverse of on 
the background, we have Q!|_ = ± ea thus Q!|_ = 
i Ea Sij , where the additional ± comes from the branches 
of in (TTOll . From the transformation d?]), the fact 

that iV® = 0 gives (c ± ea) = 0; however, condition m 
requires that (c ± ea) 7^ 0. Thus, the case with ^® 7^ 0 
is inconsistent. Indeed, although the transformation © 
is non-linearly well defined, being perturbatively non- 
invertible when 5® 7^ 0 implies that we cannot determine 
^® in terms of the background value of the old variables. 
This leaves non-uniquely determinecQ. 

Summarising, around a Minkowski background the 
only consistent branch of the matrix '/X that preserve 
Lorentz symmetry is the usual dRGT one. On the other 
hand, we have shown that there is also a ghost-free 
branch of massive gravity associated with a flat rota- 
tionally invariant background that breaks Lorentz boosts. 
Following the same logic, other consistent branches can 
be constructed choosing less symmetric backgrounds. 

Perturbative expansion. The very same results 
are obtained starting from a perturbative definition of 
the square root around a background solution of (|1]). To 
have a well defined perturbative expansion for this class 
of potentials, it is necessary that the perturbations of the 
matrix square root can be expressed in terms of the ones 
of the original matrW i.e. that the Sylvester equation 
has a unique solutioiiJ 

SVX-\/X + \/X-SVx = SX. (17) 


Actually a formal, non-standard expansion for a/X exist, where 
the degeneracy of each order in perturbations is removed by the 
next order. The potentials however feature the peculiar presence 
of square-roots of combinations of metric fluctuations in the per¬ 
turbative expansion. We leave the development of this method 
for a future work. 

^ From now on, we will leave understood the subscript |br the 
\/X when we refer to a general branch. 


It is known 11211 . and recently brought back to the at¬ 
tention p^ . Ildj . that in order to have a unique solution 
forSVX in ea, the spectrum of the eigenvalues of \/X, 
a{VX), and of —y/X, a(—y/X), should not intersect, i.e. 

cr{y/x) n a{—y/x) = 0 . ( 18 ) 

This result selects the backgrounds around which a given 
potential can be expanded according to the Sylvester the¬ 
orem. Notice that the case ^® 7^ 0 of the previous section 
exactly violates the condition (|T^ in the background. 

For instance, if we wish to have a Minkowski back¬ 
ground for g, this implies that at the background level X 
is the identity matrix. Therefore there are only two pos¬ 
sible allowed background values for y/X consistent with 

condition (fTSl) . namely y/X |±= ±14. The resulting per¬ 
turbative construction reproduces to all orders the dRGT 
potential^. Thus, dRGT is the only potential among the 
class in Q that allows a standard perturbative expansion 
around Minkowski. _ 

On the other hand, if we require for y/X only rota¬ 
tional invariance instead of the full 5'0(3,1) invariance, 
new non-trivial different branches are possibile. Taking 
as background (El), in this case X = diag(c ^,1,1,1), 
hence up to an overall sign, we can have two differ¬ 
ent branches for y/X in accordance with (El, namely 
y/X |±= diag(± c“^, 1,1,1). While the branch y/X |+ is 

the 0-order dRGT one, the branch y/X |_ will lead to 
the new one given in (1161) . 

Of course, following the same lines, we can go further 
and consider a background where also the rotational in¬ 
variance is broken producing new ghost-free potentials 
corresponding to other branches of ((4|). Again, the idea 
consists of selecting a background value for the metric 
that removes the common eigenvalues between the spec¬ 
tra of y/X and —y/X, satisfying therefore the Sylvester 
theorem (ITS)) . We will not consider these other branches 
here, however they can be potentially interesting when 
studying Bianchi type cosmological solutions in the con¬ 
text of massive gravity. 

The new potential. Let us now focus on the po¬ 
tential lArew of (HI) realised through the branch (fTHl) (let 
us set Ea = 1). Expanding around (IT^ and setting 
dfj.i' = we get that all linear terms are absent 

(e.g. g is a solution of the equations of motion) when 

oq -|- 3 Oi -|- 3 02 T <23 = 0 ! oi -l- 2 02 -l- 03 = 0 . (19) 

Notice that c represents a fiat direction and is not de¬ 
termined by background equations. This stems from the 
fact that the lapse is a Lagrange multiplier [l^. As ex¬ 
pected, we need two tunings instead of the one of the LI 
case. Expanding the potential part S'vyiew of the action 


The minus case is equivalent to the plus one, modulo a trivial 
redefinition of the coefficients an- 
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we get 


= J (^2 + 03^) (c + 1) ^ 

( 20 ) 

It should be stressed that, although from (1^01) the limit 
c —>• 1 looks regular, this limit is in fact ill defined: in¬ 
deed in the perturbative construction of y/X singularities 
are encountered in this limit. This is in agreement with 
the result from the existence of a unique solution for the 
Sylvester equation and from the ADM analysis. There¬ 
fore, the limit c —^ 1 in the y/X and in its symmetric 
polynomials do not commute. 

Given the form of the quadratic action (EUl) . one can 
check that - due to the lack of mass term for vectors 
perturbations /ipi - only 2 dof are present in the linearized 
theory [l6l - [T§| . On the other hand we know that non- 
perturbatively 5 dof are present, therefore 2 vector modes 
and a scalar one (classified according SO(3) group) are 
strongly coupled around the flat background (IT^ . The 
vanishing of is rather general [T^ and stems from the 
Lorentz invariance of the potential V. Indeed, following 
the same method used in a classic paper by Goldstone, 
Salam and Weinberg 1131 ( see for instance chapter 19.2 of 
QFT Weinberg’s book [^) to prove Goldstone theorem, 
using a multi-index notation —>• qa-, residual LI of V 

gives 


dV 


( 21 ) 


where T is a generic generator of SO{3, 1). Differentiat¬ 
ing the previous relation with respect to gA and evalu¬ 
ating the result on the solution (IT^ of the equations of 
motion, the first derivative of V vanishes and we get 


d^V 

dgAdgs g 


Tbc = 0. 


( 22 ) 


Thus the mass matrix has a zero eigenvalue for each of 
the generator of SO{3, 1) that does not annihilate g. The 
broken generators are precisely the 3 boost while rota¬ 
tions are unbroken. As a result the hgi direction repre¬ 
sents a non trivial eigenvector with zero eigenvalue of the 
mass matrix and then = 0. 


Concerning cosmological solutions, we hnd that at the 
background level this potential shares the same unfortu¬ 
nate feature of dRGT, namely it only admits open FRW 
solutions. Also in this case however, the problem can 
be overcome in the context of bigravity |2l| . It would 
be interesting to check whether the instabilities found 

in I22II and 


[23| are present when the new potential is 
use43. Also, since the vDVZ discontinuity [2^ is ab¬ 
sent for the background ca) that spontaneously breaks 
Lorentz boosts [T^, it would be useful to check whether 
new exact spherically symmetric solutions exist as the 
ones investigated in [25 1. 


Conclusions. Building on the general canonical anal¬ 
ysis of massive gravity, we have shown that new ghost 
free, Lorentz invariant massive gravity theories exist. 
They are associated with the different branches of the 
square root of X = g~^rj that enters in the construction 
of the massive gravity potential. The dRGT potential 
is the unique that produces a Lorentz invariant expan¬ 
sion around Minkowski space. However, there are other 
branches, different from dRGT, characterised by the fact 
that the expansion around flat space breaks “sponta¬ 
neously” the Lorentz invariance. Given the number of 
parameters, a flat background exists if at least an SO{2) 
of S'0(3,l) is left unbroken. We have discussed the ex¬ 
plicit case where the residual group is SO{3). For this 
case, symmetry arguments show that the vector pertur¬ 
bations (hoi) have vanishing mass. It will be of great 
interest to study the stability of perturbations for the 
cosmological solutions of this new potential. Moreover, 
since this theory is expected not to exhibit the vDVZ 
discontinuity, it will be interesting to study the features 
of its spherically symmetric solutions. Finally, it will be 
important to investigate the ghost free character of the 
new branches using the first order formalism. 
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